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Recent timing observations have revealed the presence of orbital precession due 
to spin-orbit coupling and rapid orbital decay due to dynamical tidal interaction 
in the PSR J0045-7319/B-star binary system. They can be used to put concrete 
constraints on the age, initial spin and velocity of the neutron star. 
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Introduction 

One of the fundamental questions in the studies of pulsars concerns the physical 
conditions of neutron star at birth. Of particular interest is the initial spin 
periods and velocities of pulsars, as they are related to such issues as supernova 
explosion mechanism and gravitational wave emission from core collapse. The 
PSR J0045-7319 binary (containing a 0.93.S radio pulsar and a massive B-star 
companion in an eccentric, 51 days orbit □) is unique and important in that 
it is one of the two binary pulsars discovered so far that have massive main- 
sequence star companions (The other one is PSR B1259-63). These systems 
evolve from MS-MS binaries when one of the stars explode in a supernova 
to form a neutron star. Thus the characteristics of such pulsar binaries can 
potentially be used to infer the physical conditions of neutron star formation. 
The PSR J0045-7319 system, in particular, owing to its relatively small orbit 
and "clean" environment (the mass loss from the B-star is negligible) , exhibits 
interesting dynamical orbital behaviors, which allow for concrete constraints 
on the initial spin and kick of the pulsar. 



Dynamics of PSR J0045-7319/B-star Binary 

Spin-Orbit Coupling: Based on earlier timing data, it was suggested 
that classical spin-orbit coupling is observable in the PSR J0045-7319/B-star 
binary □. This spin-orbit coupling results from the flattening of the B-star due 
to its rapid rotation: the dimensionless distortion is related to the spin rate fl s 
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by e ~ Slg/(GM C /Rc), and the resulting quadrupolc moment is Q ~ kAI c R^.e, 
where M c , i? c are the mass and radius of the B-star, k is a constant measuring 
the mass concentration inside the star. There are two effects associated with 
this spin-induced quadrupole moment: (i) The advance of periastron due to 
perturbation in the interaction potential AV ~ GM p Q/r 3 (where M p is the 
pulsar mass, r is the separation); (ii) When the B-star's spin S is misaligned 
with the orbital angular momentum L, there is an interaction torque N ~ 
(GM p Q/r 3 ) sm 9 (where 9 is the angle between S and L) between the spin 
and the orbital motion, giving rise to precessions of S and L around a fixed 
J = L + S. Both of these effects have been confirmed by recent observational. 

Rapid Orbital Decay: Recent timing data also reveal that the orbit is 
decaying!], on a timescale of P or b/-forb = —0.5 Myr (shorter than the lifetime 
of the 8.8Mq B-star and the characteristic age of the pulsar). Since mass loss 
from the B-star is negligible (as inferred from dispersion measure variation), the 
orbital decay must have a dynamical origin. It was suggested that dynamical 
tidal interaction can do jobel: Each time the pulsar passes close to the B-star, it 
excites internal oscillations (mainly g-modes) in the star, transferring orbital 
energy to the stellar oscillations. An interesting prediction of this theory is 
that in order for the energy transfer to be sufficiently large to explain the 
observed orbital decay rate, S and L must be not only misaligned but also 
more or less anti-aligned. The reason that retrograde rotation can significantly 
increase the tidal strength is the following: During a periastron passage, the 
most strongly excited modes are those (i) propagating in the same direction 
as the orbital motion, (ii) having frequencies in the inertial frame comparable 
to the "driving frequency" (equal to twice of orbital frequency at periastron) , 
and (iii) coupling strongly to the tidal potential. Since the higher-order (lower 
frequency) g-modes have smaller couping coefficients than the low-order ones, 
the trade-off between (ii) and (iii) implies that the dominant modes in energy 
transfer are those with frequencies higher than the resonant mode. If the B- 
star were nonrotating, the dominant modes would be gs-gg, and the inferred 
Poib would be two orders of magnitude too small to explain the observed 
value. However, a retrograde rotation "drags" the wave modes backwards and 
reduces the mode frequencies in the inertial frame. As a result, energy transfer 
is dominated by lower-order modes, which couple much more strongly to the 
tidal potential. At Q s — — 0.4(GM/i? 3 ) 1 / 2 (projected along the L axis), for 
example, the dominant modes are g3-gs, and the energy transfer increases 
by two orders of magnitude as compared to the nonrotating value. A recent 
analysis El of the radiative damping of g-modes indicates that an additional 
ingredient, i.e., differential rotation, is needed to compensate for the longer 
damping times of lower-order g-modes. 
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Constraints on the Initial Conditions of PSR J0045-7319 



Evidence for Supernova Kick: As mentioned before, the PSR/B-star 
binary evolves from a MS-MS binary. At this earlier stage, the B-star's spin is 
most likely to be aligned with the orbital angular momentum. The only way to 
transform this aligned configuration into the current misaligned configuration 
is that the supernova was asymmetric and gave the pulsar a kick Em, and the 
kick velocity must have nonzero components (i) in the direction out of the 
original orbital plane, and (ii) in the direction opposite to and with magnitude 
larger than the original orbital velocity. Let the total mass, semimajor axis, 
eccentricity of the system before and after the supernova be (Mj, a*, 0) and 
(Mf, a/, e/). The kick velocity is given by 

\V\ = {GMf/af) 1 ' 2 [2£ - 1 + ^ - 2(1 - e^/^V 172 costf] , (1) 

where r/ = M//M* < 1 and £ = a f / ai . With (1 + e/)" 1 < f < (1 - e,/)" 1 
and 125° < 9 < 155° (as constrained by the measurement of precessions u and 
for retrogarde rotation; tighter constraint can be obtained using the observed 
surface velocity of the B-star, but it depends on the precession phase), we get 
\V\ £ (GMf/affl 2 ~ 125 kms" 1 , where we have used the current observed 
values for a/ and e/. Orbital evolution since the supernova tends to make a/ 
and ef larger, hence decreases this lower limit. 

Age of the Binary and Initial Spin of the Pulsac: The theory of 
dynamical-tide induced orbital decay gives a scaling relationta 

Po r bOcF OTb 7/3 - 4 ^(l-e)- 6 ^), (2) 

and a similar relation for e, where v lies in the range 0.2 — 1.0, reflecting 
the uncertainty in the rotation rate. The proportional constant depends on 
the (uncertain) mode damping time. Using the observed value of P rb for the 
current system, this constant can be fixed, and the equations can be integrated 
backward in time. It was found that regardless of the uncertainties, the age of 
the binary since the supernova is less than 1.4 Myr. This is significantly smaller 
than the characteristic age (3 Myr) of the pulsar, implying that the latter is 
not a good age indicator. The most likely explanation for this discrepancy is 
that the initial spin period of the pulsar is close to its current value. Thus the 
pulsar was either formed rotating very slowly, or has suffered spin-down due 
to accretion in the first ~ 10 4 years (the Kelvin-Holmholtz time of the B-star) 
after the supernova (E. van den Heuvel, private communication). 
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